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We perform a comprehensive study of stability of a pumped atom laser in the presence of pumping, 
damping and outcoupling. We also introduce a realistic feedback scheme to improve stability by 
extracting energy from the condensate and determine its effectiveness. We find that while the 
feedback scheme is highly efficient in reducing condensate fluctuations, it usually does not alter the 
stability class of a particular set of pumping, damping and outcoupling parameters. 
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I. INTRODUCTION 

The experimental realization of Bosc-Einstein conden- 
sates (BECs) has provided a testbed for many fundamen- 
tal issues in interacting quantum systems, as well as pro- 
viding a general tool for investigating aspects of atomic 
physics such as the behaviour of weakly interacting al- 
kali gases 1]. One major application that BECs offer is 
the possibility of creating an atom laser [2| . Just as the 
optical laser revolutionised optics by offering spatial and 
temporal coherence, high spectral density and mode se- 
lectivity, the pumped atom laser offers the possibility of 
doing the same for atomic physics. 

Atom lasers have been achieved by outcoupling atoms 
from a BEC using some external means to change the 
state of a subset of the atoms in the condensate from 
a trapped to an antitrapped state El 0- This can 
produce a beam of atoms that exhibits both spatial and 
temporal coherence 0,0- As in optical lasers, a narrow 
linewidth (i.e. small momentum spread) can be attained 
by ensuring the outcoupling is weak, although this results 
in a low beam flux y| ■ In order to create a high flux, nar- 
row linewidth atomic beam, there must be competition 
between a depletable pumping mechanism and damping 
resulting in gain narrowing |7fl. 

A multimode analysis of an atom laser is possible us- 
ing semiclassical techniques such as the Gross-Pitaevskii 
equation |lj, |p| . Such a model can describe spatial varia- 
tion of the modes and determine if the laser approaches 
single-mode operation. It cannot, however, calculate the 
linewidth of the laser in this limit, since the linewidth of 
a single-mode laser is governed by the quantum statistics 
of the mode, so that a full quantum analysis is required 
[To| . This drastically increases the difficulty of the calcu- 
lation, limiting current efforts to only a few modes. 

In order for an atom laser to reach single mode op- 
eration it is necessary that it is stable, such that the 
dynamics of the condensate has a steady-state solution. 
It has previously been shown that if the pumping of the 
condensate is spatially uniform, then an atom laser will 
be unstable if the nonlinear interaction strength between 
the atoms is below some critical value This is be- 

cause increasing nonlinearity causes greater damping of 
higher frequency trap modes, leading to the promotion of 



lower trap modes, and ideally the ground state. Stability 
can also be improved if a spatially dependent pumping 
scheme is utilized, where the condensate is preferentially 
pumped in a narrow range near the center of the trap 
|l2j . As lower frequency modes have a greater overlap 
with the central portion of the trap, these modes are 
more effectively populated than higher frequency modes, 
again leading to the promotion of low-frequency modes 
and encouraging stability. 

Although interatomic interactions can lead to single 
mode operation of an atom laser, strong interactions 
will ultimately limit the linewidth, as they will cause 
phase diffusion of the lasing mode |13L fl4| . An indepen- 
dent method of improving modal stability without using 
high interaction energies is to use a feedback mechanism 
whereby energy is removed from the condensate using 
continuous knowledge of the atomic cloud's dynamics to 
tune the trap parameters. It has been shown that if 
one considers an isolated condensate in a trap, with no 
pumping, loss or outcoupling present, then provided the 
quantities (x), {x 2 ) and (IV'I 2 ) can be measured, where 
(q) = J ip*qipdx, then it is almost always possible to 
extract energy from the condensate, and thus draw it to- 
wards the lowest energy ground state [l5|. It is possible, 
however, that some specific (non-ground) state may exist 
that provides no error signal to the feedback loop causing 
the feedback process to have no effect. This occurs when 
all the moments chosen as error signals are stationary, 
i.e. do not change with time. 

The purpose of this paper is to tie these threads to- 
gether into a comprehensive analysis of atom laser sta- 
bility, incorporating loss, outcoupling, spatially depen- 
dent pumping, a depletable reservoir and feedback. We 
will consider stability across a broad range of pumping 
regimes, and include the effects of the outcoupled beam 
on the condensate as this modifies some of the stability 
conclusions in Ref. Hal. 



II. MODEL 

As discussed in the introduction, determining the sta- 
bility of an atom laser does not require a full quantum 
analysis. We will therefore use a multimode semiclassi- 
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cal model based on the Gross-Pitaevskii equation. We 
denote the condensate field by ip t (x) and the untrapped 
field by ip u (x). ^t(x) forms the lasing mode, and is 
pumped by an incoherent reservoir of atoms which has a 
density described by n(x), with the coupling between the 
two given by k(x). For calculational tractability, we re- 
strict ourselves to a one dimensional condensate, so that 
the dynamical equations for the fields can be written as 
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where m is the atomic mass, Vt is the trapping poten- 
tial, g is the acceleration due to gravity (taken to be in 
the negative x direction), Uij = ^Trh 2 aij/m is the in- 
teratomic interaction between ip t and ipj and ay is the 

s-wave scattering length between those same fields. 7. 

is the loss rate of ipi due to background gas collisions, 

(2) 

7- ' is the loss rate of tpi due to two-body inelastic colli- 

sions between particles in that state, Jtu is the loss rate 
of each field due to two-body inelastic collisions between 
particles in the other electronic state, K ut is the coupling 
rate between the trapped field and the output beam, k is 
the momentum kick due to the coupling process, r is the 
rate of density increase of the incoherent cloud of atoms 
forming the reservoir, j p is the loss rate of that cloud 
and A is the spatial diffusion coefficient. We choose the 
coupling between the reservoir and the trapped field to 
be of the form 



k(x) = KqC 
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enabling us to consider a spatially dependent pumping 
scheme. 

The pumping terms in the above equations are phe- 
nomenological, describing an irreversible pumping mech- 
anism from a reservoir which can be depleted but is re- 
plenished at a steady rate. These two features are nec- 
essary for any pumping mechanism that generates gain- 
narrowing through the competition of the gain and loss 
processes of the lasing mode. We have not included three- 
body losses, which can be important. Near a Feshbach 
resonance they may be negligible, however, and still allow 
a wide range of scattering lengths 0] . 

To model the feedback stabilization scheme we adopt 
our previous approach |15| , and assume that we can con- 
trol the condensate in some fashion in response to a set 
of continuously measured error signals. Realistic con- 
trol parameters we can use to affect the condensate will 



be the position of the trap minimum potential, the trap 
strength, and the nonlinear interaction strength. Conse- 
quently we assume it is possible to control an external 
potential of the form Vj b — a\(t)x + ci2(t)x 2 and a non- 
linear interaction strength of the form V 2 b — b(t)\%p t \ 2 ■ 
The a\(t), d2(t) and b(t) correspond to time-dependent 
control parameters that can be manipulated according 
to the measured error signals. Altering a\ and ei2 cor- 
responds to changing the position of the trap minimum 
and its curvature, while tuning b(t) corresponds to ma- 
nipulating the nonlinear interaction strength between the 
atoms in the trap. The latter can be accomplished by 
controlling magnetic field close to a Feshbach resonance 
[l7|. This is equivalent to controlling the bias magnetic 
field in a magnetic trap, or applying a constant magnetic 
field in an optical trap, and has been achieved with high 
precision [l]|. Provided the parameters ai, a% and b are 
chosen correctly, the rate of change of energy of the con- 
densate must always be non-positive. This conclusion 
assumes that the condensate is not pumped and has no 
losses. As we shall see later, the presence of such features 
can mean that this feedback scheme is not always guar- 
anteed to remove energy from the condensate, although 
in practice it normally still does a very good job. 

For our error signals we have chosen the moments (x) 
and (a; 2 ), corresponding to the first and second posi- 
tion moments of the condensate, as well as the moment 
{\4>t\ 2 ): which we have dubbed "pointiness" . 

In a real system feedback is likely to be limited due 
to finite detection speed and the ability to dynamically 
modify the potentials. As with all oscillatory systems 
controlled with feedback, when the response time of the 
feedback becomes a significant fraction of the smallest 
timescale in the dynamics of the system, the control may 
operate as positive feedback. For this reason, it is only 
safe to use control system where the dynamics of the 
relevant fluctuating moments arc within the bandwidth 
of the feedback. For most BEC systems this is not a 
problem as a control bandwidth of kiloHertz should be 
sufficient to respond to fluctuations in the system. The 
key difficulty in applying feedback in an experimental 
situation will be the destructive effects of the continuous 
detection. It has recently been demonstrated that optical 
detection cannot be arbitrarily non-destructive, and that 
for a given atomic spontaneous emission rate, increasing 
sensitivity over standard techniques requires multi-pass 
interferometry p^.l20|. 

We solved equations Q-© both with and without 
feedback terms using a pseudospectral method with a 
fourth-order Runge-Kutta time step |2l| with the XMDS 
numerical package 22] using the atomic properties and 
loss rates for 85 Rb near a Feshbach resonance, where 
the interatomic interaction can be tuned with magnetic 
fields. We used the following physically reasonable pa- 



rameters for all calculations: 7 t 
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7.0 x 10 -3 s - \ 



7t (2) = 1.7 x 10- 
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Kout = 300, k = 4.2 x 10- 4 ms _1 , A = 0.01ms- 1 , 
k = 1.0 x 10 m , with a harmonic trapping poten- 
tial of the form Vt = moj 2 x 2 /2 where ui = 50 rad s _1 
and m — 1.4095 x 10~ 25 kg. The simulation region was 
2.7 x 10 _4 m in length. As the interatomic interaction 
strengths between different Zeeman levels are unknown 
for 85 Rb we assume U u = U uu = 2U tu - 



III. STABILITY ANALYSIS WITHOUT 
FEEDBACK 

First we will consider an analysis where feedback is ab- 
sent. Previous analysis has shown that the general rule of 
stability for pumped atom lasers is that they are unsta- 
ble when the interatomic nonlinear interaction strength 
is below a certain critical value, and become more stable 
as the nonlinear interaction strength increases [llj . Fur- 
thermore, the onset of stability occurs earlier (i.e. at a 
lower nonlinear interaction strength) if the condensate is 
preferentially pumped towards the center of the trap; the 
more narrowly defined the pumping region in space, the 
more stable is the laser |T^ . 

Our more detailed analysis broadly confirms these find- 
ings. The general behaviour of the condensate is as fol- 
lows: the trap rapidly fills with atoms, reaching a cer- 
tain population governed by a balance between pumping 
and loss. The population and condensate density are ini- 
tially fluctuating, and these fluctuations will either grow 
in magnitude or become damped out over time, depend- 
ing on whether the condensate exhibits long-term sta- 
bility or not. An example of a borderline stable case is 
displayed in Figure ^ which shows the dynamics of the 
trapped (condensate) and untrapped (outcoupled) field 
densities over a period of two seconds. The condensate 
population and density rise rapidly, and then fluctuate 
about a steady state value, and then slowly stabilize. 
The density of the outcoupled field is highest near the 
outcoupling region, and decreases monotonically in the 
negative z direction, reflecting the fact that the atoms 
are accelerating under the influence of gravity. The den- 
sity fluctuations of the outcoupled field mimic those of 
the condensate. 

The dependence of stability on the interaction strength 
is also clear in our simulations, as demonstrated in Fig- 
ureEl The plots show increasing stability as the nonlinear 
interaction strength between the atoms is increased. The 
fact that the dynamics of the untrapped field closely fol- 
lows the fluctuations of the trapped field is again clear. 
Consequently we need only examine the stability of the 
condensate in order to determine the stability of the atom 
beam. 

Determining absolute stability from our simulations 
can be difficult, particularly in borderline cases where 
long timescales are required in order to determine the 
asymptotic behaviour of the system. A good technique 
for determining stability is to decompose the density fluc- 
tuations of the condensate into their Fourier components 




FIG. 1: General behaviour of the condensate and outcoupled 
fields. Shown are the condensate density over a two second 
period (top) and the outcoupled field density over the same 
period (bottom). Note that the fluctuations of the untrapped 
field follow those of the trapped field. System parameters 
were a = 1.0 X 10" 10 , a = 9.0 x 10" 6 m, r = 3.7 x 10 8 s _1 . 



and performing a modal analysis. The ground state of 
the condensate has the minimum energy allowed and is 
stable in the trap. As the energy of the condensate is 
increased, it acquires components of higher order modes, 
for example breathing modes or sloshing modes. Each 
mode will experience and gain and loss at a different 
rate, so between absolute stability and absolute instabil- 
ity exists a regime where the energy of only a subset of all 
available modes is increasing. This regime is ultimately 
unstable, but can provide a useful distinction. 

Figure [3] shows an example of the Fourier technique. 
The quantity of interest is the density of the condensate 
at the center of the trap. A time series of this quan- 
tity over a two second simulation is taken, and a Fourier 
decomposition is made of the two periods 1.0s— 1.5s and 
1.5s-2.0s. Comparison of these two results can determine 
which modes are gaining energy and which are losing en- 
ergy. In the case of Figure |3 the energy in every mode 
has decreased in the second time series, indicating that 
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FIG. 2: Increasing the nonlinear interaction leads to stabil- 
ity. Left, trapped field density; right, untrapped field density. 
Darker regions indicate higher densities. From top to bot- 
tom the scattering lengths are a = 0, a = 1.0 x 10 _10 m, 
a — 4.65 x 10 _10 m. Other parameters were a — 9.0 x 10 _6 m, 
r = 3.7 x 10 8 s~\ 
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FIG. 3: Central density of the condensate over a two second 
period (top), and the associated frequency power spectrum 
over the period 1.0s-1.5s (bottom left) and 1.5s-2.0s (bottom 
right). 



for this set of parameters the condensate is stable and is 
attracted to the ground state. 

Our previous analysis neglected the coupling between 
the trapped field and the untrapped field on the grounds 
that the effect of this coupling on the dynamics was 
usually negligible. In this approximation the dynamical 
equations have even parity, and consequently only even 
modes can be excited. The odd parity of the gravita- 
tional potential means that odd excitations can be ex- 
cited when the coupling is included. The odd modes can 
have a higher gain to loss ratio than the even modes, so 
including the outcoupling can affect the details of the be- 
haviour of the model. This is particularly evident near 
the border of stable and unstable behaviour, where the 
odd modes can grow over time resulting in instability, 
where a model without outcoupling would predict stabil- 
ity. The precise boundaries depend sensitively on the in- 
teraction of many variables, including the pumping rate. 
As we saw in the absence of spatially selective pump- 
ing |ll| , our simulations show that higher pumping rates 
lead to greater stability. A stability phase diagram is 
shown in Figure 0] showing how the boundaries of stabil- 
ity and instability are affected by altering the pumping 
rate. While the long term behaviour becomes asymptot- 
ically difficult to integrate near the stability boundaries, 
the resolution in our simulations is greater than the size 
of the features in those boundaries. This means that the 
stability boundaries are not smooth, and are very sensi- 
tive to the details of the system. 

There is also the interesting possibility that there exist 
stable condensate configurations that are not the ground 
state. For example, Figure|SJa) shows an example where 
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FIG. 4: The effect of the pumping rate on the stability of the 
laser when feedback is not present. Above the upper line in 
each plot the laser is absolutely unstable; below the lower line 
it is absolutely stable. Between the lines only certain modes 
are stable. Shown are pumping rates of r = 3.0 x 10 7 s _1 (top 
left), r = 1.0 X 10 8 s _1 m (top right), r = 6.0 x 10 8 s _1 (bottom 
left) and r = 3.0 x 10 9 s _1 (bottom right). 
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FIG. 5: An example of a set of parameters where the attractor 
for the condensate is not the ground state, but rather the 
first excited state with energy 3/2?kj per particle. The top 
figure shows the density distribution of the condensate over 
time, and the lower figure shows the energy per particle in 



the condensate, 
r = 3.7 x 10 s . 



Parameters are a = 0.cr = 10. x 10 b , and 



the trapped condensate is attracted to the first excited 
trap mode. As this represents a case where the scattering 
length is zero, the energy of the condensate has no non- 
linear component and is simply 3/2fno, the first excited 
state of the harmonic oscillator. 



IV. STABILITY ANALYSIS WITH FEEDBACK 

Stabilising the condensate can also be achieved by re- 
moving energy from it via feedback. If enough energy can 
be removed the condensate will be in the ground state 
and stable, and single-mode operation of the atom laser 
will automatically follow. We will treat the feedback as 
a set of external potentials which we can continuously 
modify based on continous measurement of various con- 
densate properties. Ideally we wish to show that if given 
a certain set of controls, it is possible to adjust them in 
such a way that the energy of the condensate will always 
be reduced. As explained in Section^ it is necessary that 
we can measure these quantities and modify them on a 
timescale shorter than the shortest relevant timescale of 
the moments we are using as error signals. 

We begin by considering the case where the conden- 
sate is isolated, so that there is no pumping, damping or 
outcoupling present. In this case condensate evolves via 



(# + Vf 6 )VM) 



(5) 



where Vfb describes the external feedback potential. Hq 
describes the evolution of the condensate when feedback 
is not present and is given by 
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where Vq is the trap potential and Uq the nonlinear in- 
teraction strength. Defining the energy of the system as 
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it is possible to show that once feedback is applied, the 
rate of change of energy is given by [l5l | 



dEo 
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As mentioned in Section [i] we choose our feedback poten- 
tial to be given by 

Vfi, = <n(t)x + a 2 {t)x 2 + b{t)\^ t \ 2 . (9) 

If the parameters a\(t), ci2(t), and b(t) are chosen to be 
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where c\, ci and C3 are positive constants and (q) ■ 
J iptQ^tdx, then the rate of change of energy becomes 
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(13) 

This is clearly non-positive and demonstrates that in this 
system the energy of the condensate can always be re- 
moved with this feedback scheme. 

In what follows, we will choose the constants to be 



ci = 2y m(muj 2 + 2cid{x 2 )j dt) 
C2 = m 2 uj 2 /h 
u\ = h 2 /mu>N, 



(14) 
(15) 
(16) 



where N is the number of particles in the condensate. 
The values of ci and C2 are chosen to ensure critical 
damping of trap oscillations in the absence of any non- 
linear interaction between atoms, and U\ is chosen such 
that it is an experimentally feasible nonlinearity that is 
efficient in removing energy from a condensate in the ab- 
sence of pumping and damping. On the occasions where 
the term proportional to C2 in l|14|) is negative and large 
enough to cause the quantity under the square root to be 
less than zero, C2 is set to zero. 

The previous analysis becomes more complex, how- 
ever, when pumping and loss terms are introduced. Now 
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when the dEo / dt is calculated one finds that there are two 
seperate contributions, one involving the feedback poten- 
tial and the other involving the pump and loss terms. 
These two contributions are decoupled; that is there is 
no term involving a combination of Vfb and pump or loss 
terms. Consequently the two can be considered inde- 
pendently. The first contribution is due to the feedback 
potential, and in the notation of Equations and © 
is given by 



dEo, 
dt 



\fb 
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ih 
2m 



The second part of i|17fl is equal to © , which is known to 
be always non-positive in our scheme. The first part of 
(|17|l represents the interaction between the outcoupling 
and the feedback, and can be positive or negative depend- 
ing on the wavefunctions of the trapped and untrapped 
fields. 

The second contribution to dE^/dt is independent of 
the feedback potential and can be broken into six parts; 

(2) 



five parts proportional to the loss terms 7^ , 7u 



7t 



(2) (2) 

7it , and j tu , and one part proportional to the pumping 
rate K,(x)n(x). The exact form of the terms is lengthy 
and will not reproduced here. The crucial point is that 
some of terms can be positive depending on the form 
of the trapped and untrapped fields. Consequently the 
feedback cannot be guaranteed to reduce the energy of 
the condensate in all circumstances, although the smaller 
the pumping, outcoupling and loss terms, the better it 
will do. 

Although the feedback scheme does not produce sta- 
ble steady state operation in all parameter regimes, it 
is still remarkably efficient at stabilising fluctuations in 
the condensate. Figure shows an example of feedback 
being applied to a condensate which is in a region of pa- 
rameter space where it is normally highly unstable. The 
introduction of feedback will often damp the condensate 
fluctuations by up to two orders of magnitude or more, 
making it useful where all that is required is an "effec- 
tive" stability requirement. That is, if one requires the 
atom laser to be in single mode for a short period of 
time, say seconds to minutes, then this feedback scheme 
can be very useful. If, on the other hand, the required 
criterion is absolute stability, then this scheme offers a 
marginal improvement. In general, the introduction of 
feedback does not greatly alter the stability class of a 
particular set of parameters. If a condensate is unsta- 
ble, such that there exist non-ground state modes whose 
energy increases with time, then applying feedback will 
drastically damp the oscillations and slow their rate of 
growth, but given time the condensate will still show the 
same level of fluctuations as before. Nevertheless, there 
are regions of parameter space which show an absolute 
improvement in stability when feedback is applied. Fig- 
ure shows an example of such a case. 
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FIG. 6: The effect of feedback. Top row: central density of 
the condensate and its energy when no feedback is present. 
Bottom row: central density of the condensate and its energy 
when feedback is turned on at t = 0.3s. Note that ground 
state energy for this system is 0.839hu> per particle after the 
particle number has reached steady state (calculated using an 
imaginary time algorithm). Parameters used are r = 3.7 x 
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FIG. 7: How adding feedback improves stability. Dashed 
line is the boundary of stability with feedback; solid line the 
boundary of stability without feedback. Pumping rate was 
r = 6.0 x 10 7 s -1 . 



V. CONCLUSION 

We have numerically simulated a pumped atom laser, 
taking into account the back action of the outcoupled 
beam as well as a variety of loss mechanisms, in order to 
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determine what factors act to stabilize the laser. As well 
as considering the effects of altering the pumping enve- 
lope, the pump rate, and the atomic scattering length, 
we introduced a feedback scheme which dynamically al- 
ters the trap parameters and atomic scattering lengths in 
order to remove energy from the condensate and reduce 
fluctuations. 

As noted in previous work, the three significant de- 
terminers of stability are the atomic scattering length, 
the pumping rate and the shape of the pumping enve- 
lope. This does not change significantly in the presence 
of the feedback scheme. Stability increases with scatter- 
ing length and pumping rate, and also increases as the 
width of the pumping region is decreased. The latter 
occurs because a narrow pump region preferentially ex- 



cites lower energy modes at the expense of higher energy 
modes. In contrast to previous work, we included the 
effect of the backaction of the outcoupled beam on the 
condensate, and demonstrated that odd modes are the 
first to become unstable. 

Introducing a feedback scheme had mixed results. It is 
certainly highly efficient at damping fluctuations in the 
condensate, leading to approximately single-mode oper- 
ation. However in many cases it does not change the sta- 
bility class of a set of parameters. That is, if a particular 
combination of pump rate, pump envelope and scatter- 
ing length is known to lead to an unstable condensate, 
where the fluctuations grow with time, applying feedback 
will drastically reduce the rate at which the fluctuations 
grow, but the system is still ultimately unstable. 
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